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Please replace the first paragraph on page 12 with the following rewritten paragraph: 



Referring to FIG. 3, uniformly distributed random conformations are generated 
allowing only rotatable bonds to vary 310. For example, 1,000 uniformly distributed random 
conformations can be generated varying only the rotatable bonds. The internal energy of each 
conformation is then minimized, again allowing only rotatable bonds to vary 320. Internal 
energy can be estimated, for example, using van der Waals potentials and dihedral angle 
term, reference: Diller, D.J. and C. L.M.J. Verlinde "A Critical Evaluation of Several Global 
Optimization Algorithms for the Purpose of Molecular Docking,'* Journal of Computational 
Chemistry, 1999, Vol. 20(16), p. 1740-1751, which is hereby incorporated herein by reference 
in its entirety. Each conformation can be minimized using, for example, a BFGS (Broyden- 
Fletcher-Goldfarb-Shanno) optimization algorithm. The algorithm is an updating formula used 
to iteratively build up an approximation of the minimum of a function f. The formula has the 
form: 



'(v/.^^-vy;.).^,..(v/^^^-vy;.)]«^« 



where u is defined as the vector 
where 

Xj is the vector of the initial position of the ligand; 
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Xj is the vector of the initial position of the ligand; 

X; ^ ^ is the updated vector of the position of the ligand; 

fi is the function to be minimized, in this case, internal energy; 

V/; is the gradient, or second derivative, of the function; 

Vfj^^ is the updated gradient of the function; 

Hi is the ith approximation to the inverse of the Hessian matrix (second derivatives of f); 

and 

Hj^.^ is the updated approximation to the inverse of the Hessian matrix. 



Please replace the last p aragraph on page 17, through line 3 on page 18: 

Each remaining match is optimized using BFGS optimization algorithm as described 
above, wherein a simple atom pairwise score is the function which is minimized 640. In one 
embodiment, the score can be modeled after the Piecewise Linear Potential (see, Gehlhaar, 
D.K., et al., "Molecular Recognition of The Inhibitor AG-1 343 By HIV-1 Protease: 
Conformationally Flexible Docking by Evolutionary Programming," Chemistry & Biology, 1995, 
Vol. 2, p. 317-324, which is hereby incorporated herein by reference in its entirety) with a 
difference being that the score used herein is preferably differentiable. For this score, all 
hydrogens are ignored, and all non-hydrogen atoms are classified into one of four categories: 

iiiei. -i ■ — ■ . 

In the Claims: 



Please amend claims 1, 4, 6, 8-11, 14, 16, 18-21, 24, 26 and 28-30 as follows: 

M -^itii — — 




(Amended) 1 . A method of docking a ligand to a protein comprising: 

performing aV)re-docking conformational search to generate multiple solution 
conformations of the ligand; 

generating a birrding site image of the protein, said binding site image 
comprising multiple hot soots; 

matching hot spots Vf the binding site image to atoms in at least one solution 
conformation of the multiple iBolution conformations of the ligand to obtain at least one 
ligand position relative to the protein in a ligand-protein complex formation; and 



